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Recently it has been shown that the non-local correlations needed for measurement-based quan- 
tum computation (MBQC) can be revealed in the ground state of the AfBeck-Kennedy-Lieb-Tasaki 
(AKLT) model involving nearest neighbor spin-3/2 interactions on a honeycomb lattice. This state 
is not singular but resides in the disordered phase of ground states of a large family of Hamiltonians 
characterized by short-range-correlated valence bond solid states. By applying local filtering and 
adaptive single-particle measurements we show that most states in the disordered phase can be 
reduced to a graph of correlated qubits that is a scalable resource for MBQC. At the transition 
between the disordered and Neel ordered phases we find a transition from universal to non-universal 
states as witnessed by the scaling of percolation in the reduced graph state. 



I. INTRODUCTION 

Quantum computers use entanglement to efficiently 
perform tasks thought to be intractable on classical com- 
puters. In one model of quantum computation, called 
measurement-based quantum computation (MBQC) [1, 
2 , the entanglement is prepared in a system of many 
particles called a resource state before computation takes 
place. Given this resource state, a quantum algorithm 
proceeds by performing adaptive, single-particle mea- 
surements, with classical processing of measurement out- 
comes. This approach is convenient for physical imple- 
mentations because single-particle operations are usually 
less error prone than entangling ones. It is also a fruit- 
ful theoretical model to investigate computationally use- 
ful phases of matter which can be studied using well- 
established methods from many-body physics. 

If a resource state is to provide a quantum speed-up, it 
must have the right kind of entanglement [3J. We call a 
resource universal [1] if we can efficiently obtain the out- 
put of an arbitrary quantum computation by performing 
single-particle measurements on it. The canonical exam- 
ple of a universal resource is the cluster state [T]. 

While it is one thing to show that a resource is univer- 
sal, for it to be viable we must also be able to prepare it 
efficiently and provide some shielding against errors. It is 
hoped that these properties can be found in natural inter- 
acting spin systems equipped with an energy gap. Find- 
ing universal resources that are natural ground states is 
interesting in its own right, because it sheds some light 
on the intrinsic computational power of natural systems. 

Unfortunately, the cluster state is not a natural ground 
state. In fact, it is impossible to have a universal resource 
of spin- 1/2 particles that is the unique ground state of 
a frustration-free Hamiltonian with only two particle in- 
teractions [SHZ]- However, this negative result does not 
hold for higher level systems. For example, a gapped 
Hamiltonian with two-body interactions involving essen- 
tially 8-dimensional systems (on a honeycomb lattice) or 
16-dimensional systems (on a square lattice) can pro- 
duce ground states that are universal resources [HI [S]. 



Moreover, the tri-cluster state [10 on spin-5/2 parti- 
cles is the ground state of a frustration-free, two-body 
Hamiltonian and is a universal resource for MBQC. How- 
ever, the Hamiltonians of both of these models lack nat- 
ural symmetries. Resources with more natural interac- 
tions based on the Affleck-Kennedy-Lieb-Tasaki (AKLT) 
state (which we define in Sec. [Tl]), have also been found. 
These models are two-body, rotationally symmetric and 
Heisenberg-like. The one-dimensional AKLT state on a 
chain [3j [11] , while not universal, can be used to imple- 
ment single qubit unitaries. Theoretical constructions 
based on the AKLT state by Cai et al. [T2] and Li et 
al. [13] were shown to be universal, the latter working 
at non-zero temperature with always-on interactions. Fi- 
nally, the two-dimensional AKLT state on a trivalent lat- 
tice is universal [TU [T3] . 

A potential difficulty with these approaches is that re- 
quiring an exact Hamiltonian to produce a ground state 
is not robust: a physical Hamiltonian will be perturbed 
from the ideal one to some degree. Hence a phase that 
is universal, rather than a specific ground state, is a 
more realistic computational resource. The computa- 
tional power of certain cluster state phases have been 
studied [TMT§] . In addition, the more natural spin-1 Hal- 
dane phase can be used as a resource to perform single 
qubit unitary operations [20 , but not arbitrary quantum 
computations. 

In this paper we investigate the computational power 
of ground states in a spin-3/2 phase of matter originally 
studied by Niggemann et al. |21] . which includes the 2D 
AKLT state. We find that a large portion of the phase 
has ground states that are universal resources, follow- 
ing similar methods to [HI [15] . The phase has several 
interesting points including a unique point where only 
projective measurements (as opposed to general POVM 
measurements) are necessary, and a transition in compu- 
tational power that coincides with the phase boundary. 

The paper is structured as follows. In section [IT] we 
describe the spin-3/2 model defined in J5JJ. The compu- 
tational power of this model is explored in section III 
by generalising me thod s used for the 2D AKLT state 
[TBI |2"2] . In section IV we highlight significant features 
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in model from the perspective of MBQC. We present our 
conclusions in section Ivl 



II. MODEL DEFINITIONS 

Consider a collection of spin-3/2 particles on a honey- 
comb lattice interacting via the Hamiltonian 

Haklt = £ > (!) 

<i,j> 

where the sum is over each pair of nearest neighbours 
and 

P S=3 _ 243 g q ,29 (o o \2 , 1 (c c \3 , _99_ j \ 
~ 1440 + 360 K°* °j) t 90W' + 1152 ' \ Z > 

projects nearest neighbours i and j onto the seven di- 
mensional subspace of total spin Stot = 3. We will call 
this model the 2D AKLT model after the authors Affleck, 
Kennedy, Lieb and Tasaki who originally proposed it [23] . 
The AKLT model can be thought of as a deformation of 
the Heisenberg model H — Si • Sj that preserves full 
rotational symmetry. Note, however, that unlike the ID 
case, the AKLT model and the Heisenberg model are 
not in the same phase: the Heisenberg model has a Neel 
ordered ground state, while the AKLT model does not. 
Thus the AKLT model is said to be in a disordered phase. 
The absence of Neel order makes it a more realistic model 
for certain systems, e.g. Bi3Mn40i2, which is a spin-3/2 
antifcrromagnet on a honeycomb lattice without Neel or- 
der [23]. The ground state of the AKLT model \^aklt), 
which we will call the 2D AKLT state, is a valence-bond 
solid, or projected entangled pair state (PEPS). Details 
of the PEPS construction of ground states are included 
in Appendix |A"| 

Niggemann et al. [21] studied a 5-parameter deforma- 
tion of the 2D AKLT Hamiltonian which is frustration 
free and whose ground state is a one-parameter deforma- 
tion of the AKLT PEPS (see Appendix[A} . This Hamilto- 
nian is still two-body nearest neighbour with summands 
that preserve two Z 2 symmetries, parity and spin flip, 
however it breaks full rotational symmetry to a U(l) 
symmetry (arbitrary rotations about the z-axis). The 
deformed even parity Hamiltonian is 

H(a) = l D ( a )i ® D ( a )j] M°) L°( a )* ® D ( a )i\ > 

<i,i> 

(3) 

where 

3 

hi,j{a) = £ \"i\\Stot = 3,m)(S tot = 3,m| , (4) 

m— — 3 

Here D{a) — diag(V3/a, 1, 1, v3/a) in the S z basis and 
the continuous free parameters satisfy: Ao, Ai, A2, A3 > 
and a can be positive or negative. In this work we focus 
on the regime where a is strictly positive so that D(a) 



is a bounded positive operator but our protocol works 
just as well for a strictly negative with the replacement 
D{a) — > D(\a\). Importantly, the ground state of H(a) 
depends only on a and the ground state of H(y/3) is 
\iPaklt}- The fully rotationally invariant interaction, 
Haklt in Eq. 0, corresponds to the choice of param- 
eters a = v3 and A m = Vim. With periodic boundary 
conditions or with open boundaries and Heisenberg in- 
teractions between spin- 1/2 particles and the edges, the 
ground state \ip( a )) is unique for a > and can be ob- 
tained simply by applying the inverse deformation to the 
2D AKLT state 

|V(o)> cx (Diar^HAKLT) ■ (5) 

Using Monte Carlo sampling, Niggemann et. al [3T] found 
the ground states had exponentially decaying correlation 
functions below a critical value of a 2 = 6.46, while were 
Neel ordered above this value. Thus, Hamiltonians in the 
< a 2 < 6.46 region are conjectured to be gapped, while 
Hamiltonians in the a 2 > 6.46 region are gapless [25 . 

We will refer to the appearance of Neel order at a 2 = 
6.46 as the phase transition in this model. We will label 
the region a 2 < 6.46 as the AKLT phase, and the re- 
gion a 2 > 6.46 as the Neel ordered phase. Note that the 
area law for entanglement holds across this phase transi- 
tion (PEPS dimension is constant), a property that can 
only occur in PEPS on graphs of dimension greater than 
one |2Jj]. We also note that Schuch et al. [57] have stud- 
ied classes of PEPS related by this type of symmetry- 
preserving deformation. 

III. MBQC USING GROUND STATES IN THE 
AKLT PHASE 

In this section we will look at how ground states in the 
AKLT phase (as defined above) can be used for MBQC. 
To do this we generalise the existing method used at the 
AKLT point [141 115) , which we will briefly review. 

A. Protocol at AKLT point 

The 2D AKLT state has been shown to be a uni- 
versal resource for measurement-based quantum com- 
putation [141 115] . We will summarize the procedure 
for measurement-based quantum computing on the 2D 
AKLT state by breaking it into two stages: reducing to 
a stochastic graph state, then using this graph state for 
computation. 

1. Reduction to a stochastic graph state 

The first stage relies on the principle of quantum state 
reduction [28 , where a resource is shown to be universal 
by proving that it can be converted into a known univer- 
sal resource efficiently by single-particle measurement. A 
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FIG. 1: Illustration of how the graph state is encoded on 
the post-filtering AKLT state. In (a) we illustrate a small 
2D AKLT state on a trivalent lattice, where each node corre- 
sponds to a spin-3/2 particle, and is labelled according to the 
filter outcome obtained. In (b), we illustrate the graph ob- 
tained by processing (a) such that nodes represent domains of 
like outcomes, and edges are bonds between domains. In (c) 
we have illustrated the encoded graph state, obtained from 
(b) by deleting edges modulo 2. Each node represents an 
encoded qubit, and each edge a graph state edge. 



three-outcome filtering measurement is performed on ev- 
ery particle. Define \m)b to be the spin-3/2 state satisfy- 
ing Sb\m)b = to | to) & where Sb is the spin-3/2 component 
along the b axis, b £ {x, y, z}, and to G {§, |, — \, — §}. 
The measurement operators for the initial filtering are 
chosen to be {F x ,F y ,F z } where 



F h = 



(!§>*< 



("ID 



(6) 



These operators satisfy the completion relation 



F£F X + F^F y + F\F Z = I and thus form a valid 



set of measurement operators, i.e., {E x ,E y ,E z } := 
{F^F x ,F^F y , FlF z } is a POVM :29]. The measurement, 
applied globally, projects each spin-3/2 system onto a 
two dimensional, or qubit, subspace. We label each par- 
ticle either X, Y, or Z according to the outcome of this 
measurement. The resulting collection of spin-3/2 par- 
ticles encodes a graph state, which can be proven using 
the stabilizer formalism [14] or by using a tensor network 
description [L3] , 

The graph state is encoded as follows (we have illus- 
trated the encoding in Fig. [7J). A domain is defined as 
a connected set of particles with the same label. Each 
domain encodes a single qubit in the graph state. An 
edge exists between two encoded qubits if an odd num- 
ber of bonds (in the original honeycomb lattice) connect 
the corresponding domains. 

We remark that the reduction of 2D AKLT state via 
a three-outcome POVM to a stochastic graph state is 
similar to the reduction of the tri-cluster state to a graph 
state described in [28] . however the graphs produced in 
the latter are deterministic. 



on a square lattice, which is itself a universal resource 
|28j . Alternatively, in Ref. [IS] , 'backbone' paths are 
identified through the graph state along which corre- 
lation space qubits can propagate and interact, again 
enabling universal quantum computation. Essentially, 
both approaches use a stochastic graph state as a re- 
source. Whether this is possible depends on the stochas- 
tic graph states having certain desirable properties. We 
will show how the same approach can be applied to de- 
formed AKLT states. 



B. Generalized reduction scheme 

Here we generalize the above method to show how de- 
formed 2D AKLT states can be reduced to stochastic 
graph states using a modified version of the {F x , F y ,F z } 
measurement. For a > 1 (we will consider the a < 1 case 
in section IV 3 1 we define three measurement operators 



F x {a) 



F y (a) 



Fz{a) 




D(a)^F x D(a) , 



D(a)^F y D(a) , 



D{a)^F z D{c 



(7) 



Numerical prefactors are included to ensure that 
F x {a)^F x {a)+F y {a^F y {a)+F z {a)^F z {a) = I. The mea- 
surement operators F x (a), F y (a), and F z (a), like F x , F y , 
F z , are projections onto qubit subspaces, up to a constant 
factor. 

The reduction procedure involves performing this mea- 
surement on every particle of the deformed ground state 
(D(o,)~ 1 )® n \'iPaklt)- The resulting state after subject- 
ing every particle of the deformed AKLT state \if)(a)) to a 
POVM measurement {F x (a), F y (a), F z (a)} is equivalent 
to the state obtained by measuring the undeformed state 
\iPaklt) using the POVM {F x , F y , F z } and getting the 
same outcomes, up to local unitaries. 

Thus we can apply the existing methods in section 
|III A2| to resulting stochastic graph states away from the 
AKLT point. However the success of these methods de- 
pends on the stochastic graph states having certain prop- 
erties. The statistics that determine these properties are 
dependent on the value of a, as we will explain in the 
following section. 



2. Using stochastic graphs state as resources 

In the second stage of this method, the stochastic 
graph state is used for MBQC. Following Ref. [14], arbi- 
trary quantum computations may be performed by first 
converting the post-filter graph state into a cluster state 



C. Statistical model 

Because each particle is measured with a three- 
outcome POVM, the total number of possible outcomes 
is 3 N where N is the number of spin-3/2 particles. Some 
of these outcomes correspond to computationally useful 
graph states (e.g. if every domain had size one), while 
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FIG. 2: A typical reduction outcome at the AKLT point. In 
(a) each node corresponds to a spin-3/2 particle, and edges are 
drawn between nearest neighbours. The nodes are coloured 
according to which outcome was obtained: Z outcomes are 
cyan, X outcomes are magenta and Y outcomes are yellow. 
In (b) the resulting graph state is drawn. Each node corre- 
sponds to a qubit (a domain of like measurement outcomes), 
and edges correspond to graph state edges. Periodic bound- 
ary conditions are imposed and, for clarity, some vertices and 
edges have not been displayed. The graph has many crossings, 
making it useful for MBQC. 



some will not (e.g. if every measurement outcome was 
Z). Let a = <ti, . . . , <7jv be a sequence of filter outcomes 
where cr, is the filter outcome on spin i and is either X, 
Y or Z. At the AKLT point it was shown in [14] that 
the probability of obtaining a particular a is 



p(a) = I 2 I^)I-I S (-)I , 



(8) 



where |V(ct)| is the number of domains for a given out- 
come, |-E(<t)| is the number of intcr-domain bonds be- 
fore deleting edges in the reduction to a graph state, 
and Z = 53 , 2l y ( <T >\~\ E (" )l is a normalisation factor. 
A typical filter outcome, sampled from this distribution, 
is shown in Fig. [2j 

In Appendix[B]we explain how to use Eq. Q to compare 
the norms (hence probabilities) of the post-filter states at 
a v3 to those at a = y3- The probability of obtaining 
a particular filter outcome a with deformation a is 



p(a, a) 



Z{a) 



2 \V{cr)\-\E(a)\ 



where |V(c)| and |£7(<r)| are as above, N z (a) is the total 
number of Z filter outcomes. These statistics are equiva- 
lent to a Potts-like spin model in the canonical ensemble 



p(a,a) = 



1 



Z(a) 



,-«B( a )-VW-B(a)JT,W) 



where the E{a) term is the Potts Hamiltonian 30 , V(a) 
is a non-local cluster counting term similar to the random 
cluster model [3TJGI2], B{a)N z {a) is an external field term 
with strength B(a) = log 2 (a 2 — 1) — 1, and the inverse 
temperature (3 — log e 2 is constant. This shows that 
varying a to deform the AKLT model is like varying an 
external magnetic field in terms of the statistics of the 
filter outcomes. 



D. Identifying computationally powerful ground 
states 

Here we will show that, beyond the AKLT point at 
a = \/3j there is a range of a values that have univer- 
sal ground states. The reduction process in section [Til B| 
produces stochastic graph states with statistics given by 
Eq. ([9]) . For some filter outcomes it is possible to convert 
the stochastic graph state to a cluster state on a honey- 
comb lattice, which is itself a universal resource [33] . A 
ground state at a given value of a is universal if we can 
reduce it to a honeycomb cluster state efficiently, i.e., if 
we can produce honeycomb cluster states of size N from 
a ground state with poly(N) particles in poly(N) time. 
There are two conditions that will ensure this is possible 
QJIH]: 

1. The maximum domain size scales no faster than 
logarithmically with the lattice size; 

2. The probability of the stochastic graph state hav- 
ing a crossing (a path of edges connecting opposite 
boundaries of the graph) tends to one in the limit 
of large N. 

Condition 1 ensures that producing graph states with an 
arbitrary number of qubits is possible. It also rules out 
the possibility of an macroscopic domain, which would 
produce star-shaped graphs states (see Fig. [3J: for an ex- 
ample that is not universal for MBQC). If condition 1 
is satisfied then condition 2 will imply the existence of a 
extensive number of crossings in both lattice dimensions, 
which guarantees the existence of a honeycomb subgraph 
[To] , and hence the universality of the state. 

We performed Monte Carlo sampling over the distri- 
bution (|9) to determine which values of a correspond 
to ground states satisfying these two conditions. We 
performed simulations on lattices of varying size up to 
120 x 120 spins. Details of the numerical methods used 
are provided in Appendix [Cj Samples of resulting graph 
states are displayed in Fig. [3] We found that maxi- 
mum domain sizes scale logarithmically in the region 
(9) 1 < a 2 < 6.46 while a macroscopic (of extensive size) 
domain appears at a 2 > 6.46. The scaling of maximum 
domain size at selected values of a is presented in Fig. [4] 
and the probability of obtaining a macroscopic domain 
as a function of a is presented in Fig. [5] for various lattice 
sizes. 

To assess condition 2, we directly checked for the exis- 
tence of a crossing of the resulting stochastic graph state. 
We found numerically that the resulting graphs have 
crossings with probability one for lattices of size 20 x 20 
and larger (up to 120 x 120, the limit of our simulations), 
for all values of a. Our conclusion is that there is compu- 
tationally powerful region that ends at a 2 = 6.46 ± 0.05, 
the upper limit corresponding to the boundary between 
the AKLT phase and the Neel phase. 

Due to the presence of a macroscopic domain of fil- 
ter outcomes, our method for MBQC fails for a 2 > 6.46, 
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FIG. 3: The largest connected component of graphs sampled 
from Eq. ^ at (a) a 2 = 1, (b) a 2 = 5.70 and (c) a 2 = 6.96 on 
a 20 x 20 lattice. We use the same coloring as Fig. [2] and, for 
clarity, some edges have not been displayed. Graphs (a) and 
(b) are in the computationally useful region, while (c) is not. 
In (a) there are no Z outcomes and crossings are more sparse 
than at the AKLT point. In (b) there are large Z domains, 
which appear as vertices of high degree. In (c) we are in 
the supercritical region and there is a Z domain spanning the 
lattice. This results in a graph that is not universal for MBQC 
(it has a tree structure and cannot be efficiently reduced to 
graph state on a honeycomb lattice). 
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FIG. 5: Probability of a domain that spans the lattice vs a for 
various lattice sizes close to the critical point. The probability 
tends to a step function as N — > oo with a discontinuity at 
a = 6.46. This shows there is an macroscopic domain above 
this point, thus graph states produced above this point are 
not expected to be computationally useful. 



however we haven't ruled out the possibility that ground 
states in this region can be used as computational re- 
sources using another method. However, the universality 
of such ground states is unlikely as ground states above 
the critical point of a 2 = 6.46 are Neel ordered. While 
states with long range order are usually expected to not 
be universal for MBQC, we warn that exceptions have 
been found (34]. 



IV. EXPLORING THE PHASE 
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FIG. 4: Maximum domain size vs. number of spins in the 
ground state for selected values of a. A straight line is fitted 
to the a 2 = 6.46 data points. For values of a 2 < 6.46 domain 
sizes scale logarithmically. 



In this section we will highlight significant features of 
the model characterised by particular values of a. 



At a 2 — 3 we have the AKLT state. This point is 
optimal in the sense that it produces graph states with 
the most qubits. In contrast, as a 2 — > oo the inverse 
deformation D{a)~ 1 tends towards a projection onto the 
space spanned by |±|) z resulting in a GHZ ground state 

l/y2(|tm ■■•> + (-1)^ lltlt ••■)) where |t) = |f)zand 
|t) = | — |) z . Any measurement sequence on this state 
can be simulated efficiently on a classical computer. 



Note that F z (l) = 
surement in Eq. Q 



•■ and therefore the filtering mea- 
at a 2 = 1 has only two outcomes, 



F x (l) and F y (l). We define the orthonormal basis 



10) 


^0 




Ml), + 1-1),) 


11) 


-Id 


2)2 ~ \~ i) z ' 


HD.-l-i)*) 


|2) 


-id 






|3) 


-id 


2)2 * 1 2 ) ^ 


-|5)* + i|-D- 



Then we can write F X {1) = |0)(0| + |1)(1| and F y (l) = 
|2)(2| + |3)(3|, which are projections onto orthogonal 
spaces. Hence the a 2 = 1 ground state is special in that 
it requires only projective measurements to be universal 
for MBQC. 

3. a 2 < 1 

The filtering measurement in Eq. ([7| is not well-defined 
for a 2 < 1 . Here we will provide a casual analysis of how 
states within region may be useful. For a 2 < 1 we define 
a new measurement with the operators 

{aF x (a),aF y (a),E(a)} (12) 

where E{a) := diag(0, \/i - a 2 ,y/l - a 2 ,0). The F x (a) 
and F y (a) outcomes produce graph state qubits as be- 
fore, however E{a) outcomes must be treated separately. 
When a 2 is very slightly less than 1, the state will be like 
the a 2 = 1 state except for a few isolated E(a) sites. At 
an E{a) site we can measure surrounding X and Y qubits 
in a disentangling basis (corresponding to a Z cluster 
state measurement), effectively disentangling E(a) sites 
from the others. However, as we decrease a towards zero, 
the number of E(a) outcomes increases, and eventually 
we cannot cut them out of the lattice without adversely 
affecting the connectivity of the graph. Hence we predict 
a critical value of < a < 1 below which this measure- 
ment produces states that are not universal for MBQC. 
We leave a detailed analysis of the < a < 1 region to 
future investigation. 

V. CONCLUSION 

In summary, we have studied the computational power 
of a spin-3/2 AKLT phase that preserves U(l) and I2 
symmetries |21j . By mapping measurement outcomes 
to a classical spin model we identified three regions: a 
region with ground states that are universal resources 
(1 < a 2 < 6.46), a region that is unlikely to be computa- 
tionally powerful (a 2 > 6.46), and a region that we can- 
not say much about (0 < a 2 < 1). Significant points in- 
clude the 2D AKLT state (a 2 = 3), a state which requires 
only projective measurements (a 2 = 1), a GHZ state 
(a 2 = 00) and the phase transition (a 2 = 6.46) which cor- 
responds to a transition in computational power. While 
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FIG. 6: Illustration of PEPS ground states. Black dots rep- 
resent virtual spin-1/2 particles, lines connecting them repre- 
sent singlet bonds and large circles are locations of physical 
particles. 

it is an open question whether or not this quantum com- 
putational phase is gapped, it is known that the ground 
state for < a < 00 with periodic boundaries (or open 
boundaries with Heisenberg interactions with qubits on 
the boundaries) is unique [3T]. Any size dependent gap 
in the disordered phase would be expected to scale at 
worst as an inverse polynomial in system size N. 

A practical limitation of the method is that it depends 
on precise knowledge of the parameter a. Performing the 
procedure with an assumed value of a that differs from 
that of the actual ground state will yield a resource state 
that differs from the cluster state. The effect will be 
that X and Y outcomes cause errors in the correlation 
space in which the computation takes place (Z outcomes, 
however, are error free). It is not even clear that these er- 
rors can be corrected using standard techniques, as they 
may not correspond to linear completely-positive trace- 
preserving maps on the correlation space [35 . Whether 
there exists a method that is independent of the exact 
value of the deformation, analogous to [20], remains to 
be seen. Another question is if other deformations to 
the 2D AKLT model (e.g. ones that preserve full rota- 
tional symmetry) yield computationally powerful ground 
states. 
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Appendix A: Ground states as PEPS 

The ground states in Eq. ^ can be written as PEPS. 
We place a singlet state \ip~) = l/%/2(|01) - 1 10) ) on 
each edge of the honeycomb lattice, where |0) and |1) are 
virtual spin-1/2 states. This places three virtual spin- 
1 1 2 particles at each vertex, where a vertex corresponds 
to the location of a single physical spin-3/2 particle, as 
illustrated in Fig. [6] We obtain the physical ground 
state by applying the map D(a)T to the three spin-1/2 
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state label 


state on A sites 


state on B sites 


l*t> 


\m = 3/2) 


|m = -3/2) 


\zA) 


-\m = 1/2) 


\m = -1/2) 


\zV) 


\m = -1/2) 


|m = 1/2) 


\zl) 


-|m = -3/2) 


|m = 3/2) 



TABLE I: New basis state labels for convenience. The hon- 
eycomb lattices is bipartitioned into A and B sites, where A 
sites have a bond pointing down, and B sites have a bond 
pointing up. 



where the second term is the probability ra- 
tio at the AKLT point, shown in [22] to be 

2 \V(*)\-\E{a)\-\V(v')\ + \E(*-)\_ Thus we haye 



p{p, a ) 

p{a',a) 



N*{<t)-N z {<t') 



x 2 \V(a)\-\E(a)\~\V{<y')\ + \E(a')\ 



particles at each site where T is the projection onto spin- 
3/2. Hence the ground state can be written as 



which is equivalent to Eq. 



|V(a)> cx (9)p(a)T) t 



,|VT 



(Al) 



e£E 



which means that singlets are placed on every edge of the 
honeycomb lattice E and the projections D(a)T map the 
three virtual spin-1/2 particles at each vertex to physical 
spin-3/2 particles. 

To simplify the PEPS tensors, we define a new spin-3/2 
basis in Table |U This gives the ground states the defining 
three-index tensors 



A[zt] = |0)„/d|0)i|0) r , 
A[zA] = 1^(11)^10)^0),. 

+ |0)„ /d |l)/|0) r 
A[zV] = l/a(|l) u/d |l),|0) r 

+ |0) u/d |l),|l)r 

A[zi] = |l) u/d |l),|l) r . 



}u/d 



lu/d 



|0)«|Q)r), 



|0>z|l>r 



(A2) 

(A3) 

(A4) 
(A5) 



Appendix C: Monte Carlo sampling 



We sampled the distribution in Eq. ^ using the 
Metropolis-Hastings algorithm with single-spin flip dy- 
namics, as was done by Wei et al. P3]. We used essen- 
tially the same procedure as [M], however some changes 
were made to work with values of a 2 7^ 3. For one, we 
used ([9]) to obtain a generalised a-dependent Metropolis 
ratio, 



a? - 1 



N z (a')~N z (a) 



Appendix B: Distribution of measurement outcomes 



2 \V(a')\-\E(a')\-\V(a)\ + \E(a)\ 



We obtain the probability distribution in Eq. (J9j) by 
calculating the ratio 



P(cr,a) = yj(a)\{F a (a)}\i>(a)) 
V{<J', a) (^(o)|{^(o)}|^(a)) ' 



(Bl) 



where a and a' are two filter outcomes, and {F a (a)} = 
(a)F ai (a)®- ■ -2>F^ N (a)F aN (a). The a-dependence of 
the probability ratio is contained in the numerical pref- 
actorsofEq. ([7]), and the norms of D(a)\ ± %) x ,y,z- Using 
this we can rewrite Eq. (Bl), with the a-dependence as 
a separate factor 



p{v, a) 
p(a',a) 



a 2 -l 



Mv^)|{F g (v^)}|V>(>/3)) 

(^(V5)|{F^(V3)}|V(V3)) ' 



(B2) 
(B3) 



where a is a filter configuration in the Markov chain, 
and a' is the proposed next filter configuration (obtained 
by flipping a single spin in a). We also generalised the 
starting filter configuration to depend on a, to reduce 
burn-in time. This initial configuration was obtained by 
assigning a label (X,Y or Z) independently to each spin 
with probabilities 



Pz = 



Py = (! -P*)/2, 



(CI) 
(C2) 



where pb is the probability of assigning the label b. This 
is the probability distribution obtained by neglecting cor- 
relations between filter outcomes (the 2\ y ( cr ^ + \ E (' 7 ^ term 
in Eq. 



8 



[1] R. Raussendorf and H. J. Briegel, Phys. Rev. Lett. 86, 
5188 (2001). 

[2] H. J. Briegel, D. E. Browne, W. Diir, R. Raussendorf, 

and M. Van den Nest, Nat. Phys. 5, 19 (2009). 
[3] D. Gross, S. T. Flammia, and J. Eisert, Phys. Rev. Lett. 

102, 190501 (2009). 
[4] D. Gross, J. Eisert, N. Schuch, and D. Perez-Garcia, 

Phys. Rev. A 76, 052315 (2007). 
[5] M. Van den Nest, K. Luttmer, W. Diir, and H. J. Briegel, 

Phys. Rev. A 77, 012301 (2008). 
[6] J. Chen, X. Chen, R. Duan, Z. Ji, and B. Zeng, Phys. 

Rev. A 83, 050301 (2011). 
[7] M. A. Nielsen, Rep. Math. Phys. 57, 147 (2006). 
[8] S. D. Bartlett and T. Rudolph, Phys. Rev. A 74, 040302 

(2006). 

[9] T. Griffin and S. D. Bartlett, Phys. Rev. A 78, 062306 
(2008). 

[10] X. Chen, B. Zeng, Z. Gu, B. Yoshida, and I. L. Chuang, 

Phys. Rev. Lett. 102, 220501 (2009). 
[11] G. K. Brennen and A. Miyake, Phys. Rev. Lett. 101, 

010502 (2008). 

[12] J. Cai, A. Miyake, W. Diir, and H. J. Briegel, Phys. Rev. 

A 82, 052309 (2010). 
[13] Y. Li, D. E. Browne, L. C. Kwek, R. Raussendorf, and 

T. Wei, Phys. Rev. Lett. 107, 060501 (2011). 
[14] T. Wei, I. Affleck, and R. Raussendorf, Phys. Rev. Lett. 

106, 070501 (2011). 
[15] A. Miyake, Ann. Phys. 326, 1656 (2011). 
[16] D. E. Browne, M. B. Elliott, S. T. Flammia, S. T. Merkel, 

A. Miyake, and A. J. Short, New J. Phys. 10, 023010 

(2008). 

[17] S. D. Barrett, S. D. Bartlett, A. C. Doherty, D. Jennings, 

and T. Rudolph, Phys. Rev. A 80, 062328 (2009). 
[18] A. C. Doherty and S. D. Bartlett, Phys. Rev. Lett. 103, 



020506 (2009). 

[19] S. O. Skr0vseth and S. D. Bartlett, Phys. Rev. A 80, 
022316 (2009). 

[20] S. D. Bartlett, G. K. Brennen, A. Miyake, and J. M. 

Renes, Phys. Rev. Lett. 105, 110502 (2010). 
[21] H. Niggemann, A. Kliimper, and J. Zittartz, Z. Phys. B 

104, 103-110 (1997). 

[22] T. Wei, I. Affleck, and R. Raussendorf, arXiv: 1009.2840 
(2010). 

[23] I. Affleck, T. Kennedy, E. H. Lieb, and H. Tasaki, Comm. 

Math. Phys. 115, 477 (1988). 
[24] R. Ganesh, D. N. Sheng, Y. Kim, and A. Paramekanti, 

Phys. Rev. B 83, 144414 (2011). 
[25] B. Nachtergaele and R. Sims, Comm. Math. Phys. 265, 

119 (2006). 

[26] F. Verstraete, M. M. Wolf, D. Perez-Garcia, and J. I. 

Cirac, Phys. Rev. Lett. 96, 220601 (2006). 
[27] N. Schuch, D. Perez-Garcia, and I. Cirac, Phys. Rev. B 

84, 165139 (2011). 
[28] X. Chen, R. Duan, Z. Ji, and B. Zeng, Phys. Rev. Lett. 

105, 020502 (2010). 

[29] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
2004). 

[30] F. Y. Wu, Rev. Mod. Phys. 54, 235 (1982). 
[31] C. M. Fortuin, Physica 58, 393 (1972). 
[32] R. G. Edwards and A. D. Sokal, Phys. Rev. D 38, 2009 
(1988). 

[33] M. Van den Nest, W. Diir, A. Miyake, and H. J. Briegel, 

New J. Phys. 9, 204 (2007). 
[34] D. Gross and J. Eisert, Phys. Rev. Lett. 98, 220503 

(2007). 

[35] T. Morimae and K. Fuji, arXiv:1110.4182vl (2011). 



